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Abstract 

We consider early time cosmic evolution in string gas cosmology dominated by open strings at- 
tached to D-branes. After reviewing statistical properties of open strings in D-brane backgrounds, 
we use dilaton-gravity equations to determine the string frame fields. Although, there are distinc- 
tions in the Hagedorn regime thermodynamics and dilaton coupling as compared to closed strings, 
it seems difficult to avoid Jeans instability and assume thermal equilibrium simultaneously, which 
is already a known problem for closed strings. We also examine characteristics of a possible sub- 
sequent large radius regime in this setup. 



*ayse . aslanarginOboun . edu.tr 
^^ali .kayaSboun. edu.tr 



I. INTRODUCTION 



As a candidate for a unified theory of all interactions and quantum gravity, string theory 
still lacks any direct or indirect observational support. One hopes that cosmology may alter 
this situation and offer a testing ground for string theory or any other quantum theory of 
gravity in not so distant future. Therefore, it is important to study possible cosmological 
implications of string theory. String gas cosmology [1] is one possible approach which takes 
into account not only the massless low energy fields but all stringy excitations like winding 
modes, which are expected to play important roles both in early and late time cosmologies 
(see [2, 3] for recent reviews). 

Except in special situations usually preserving some supersymmetry, we do not have an 
understanding of string theory in the strong coupling regime. For cosmologically relevant 
time dependent backgrounds this lack of knowledge prevents proper understanding of the key 
issues like the resolution of big-bang singularity. However, this does not stop one to study 
the theory in the weak coupling limit. Indeed, string theory has a very rich structure even 
when the interactions are weak and this motivates the study of toy cosmological models. 
For instance, it is known that for strings at very high energies there exists a "limiting" 
Hagedorn temperature and one wonders possible implications of this feature in cosmology. 
On the other hand, studies in the string/brane gas cosmology reveal that stringy excitations 
offer mechanisms for late time stabilization of extra dimensions (see e.g. [4-12]), shape 
moduli (see e.g. [14-17]) and dilaton (see e.g. [12, 18-20]). 

Since the thermal partition function does not converge for temperatures above the Hage- 
dorn temperature, thermodynamics of strings should be studied in the microcanonical en- 
semble. There are subtleties in taking the thermodynamical limit in this formalism [21- 
23], but assuming a totally compact space, the microcanonical approach offers a suitable 
framework in studying string cosmology near "big-bang". To be specific, one can use basic 
thermodynamical relations to deduce the right hand side of the Einstein's equations and 
study the cosmological evolution. In string theory, at least in the weak coupling limit, all 
the necessary thermodynamical variables are already determined and ready to use for cos- 
mology. For closed strings, this approach was used in [24-31] to study different aspects of 
string gas cosmology. In [30], however, an important obstacle regarding the assumption of 
thermal equilibrium in early time string gas cosmology is pointed out. Namely, it is noticed 
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that the interaction rates for closed strings turn out to be small compared to expansion rate 
and thus strings cannot keep thermal equilibrium, which is an important barrier in such 
scenarios. 

Our aim in this paper is to study an early time cosmology dominated by open strings 
attached to D-branes. In the absence of any physical, mathematical or philosophical con- 
straints about initial conditions in the big-bang, it is natural to expect the existence of 
D-branes even at early times. Of course, according to Brandenberger-Vafa (BV) mechanism 
[1] higher dimensional D-branes are expected to annihilate, but such arguments rely on the 
assumption of thermal equilibrium which are subject to discussion as noted above. Indeed, 
for D-branes the interaction rates are studied in [32], which indicates difficulties for BV 
mechanism applied to D-brancs. 

In the presence of D-branes, it is known that open strings dominate the system in thermal 
equilibrium. Physically, this can be understood as a result of D-branes chopping up the 
closed strings [33-35]. In the micro-canonical approach, the statistical properties of open 
strings on D-brane backgrounds is determined in [35] and there are key differences compared 
to closed strings. Besides, the coupling of open strings to dilaton alters gravity field equations 
in a non-trivial way. Thus, one expects a difi^crent cosmological evolution in the presence of 
D-branes. Our aim is to check whether these differences may help to allow the assumption 
of thermal equilibrium for open strings in a simple setting. It turns out that this problem 
also persists for open strings in the weak coupling limit and thus studies of string gas 
cosmology with D-branes also require understanding of non-equilibrium thermodynamics or 
strong coupling dynamics or both. 

The plan of the paper is as follows. In the next section, following [35] , we review thermo- 
dynamics of open strings on parallel D-branes which are uniformly distributed along some 
compact directions. In section III, we use this thermodynamical information in dilaton- 
gravity field equations to determine the early time cosmological evolution. We check the 
consistency of the background and examine the reliability of the thermal equilibrium. In 
this section, we also study characteristics of a possible subsequent large radius regime. We 
conclude with brief remarks and possible future directions in section IV. 
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II. REVIEW: THE ENTROPY OF OPEN STRINGS IN D-BRANE BACK- 
GROUNDS 



In this section, we review the calculation of the entropy of open strings attached to 
D-branes in a totally compact toroidal space at temperatures close to the Hagedorn temper- 
ature and in the weak string coupling limit. The main results derived from first principles 
can be found in [35] and here wc simply repeat the derivations in a slightly different way 
for the cases we are interested in (see also [33, 34] for an approach based on Boltzmann 
equations). We assume that there are n parallel dN dimensional D-branes in the system 
which are distributed uniformly in do dimensional transverse directions, where d^ + do — 9 
(see figure 1). Actually, in the absence of any symmetry one expects to have zero net RR 
charge and thus an equal number of D-branes and anti-D-branes in the system. In the 
next section, we consider this possibility, but for now we simply assume only the existence 
of D-branes. We use Rn and Rd to denote the corresponding radii for the world- volume 
and transverse directions. The end points of open strings are imposed to obey Neumann 
boundary conditions along the world- volume directions and Dirichlet boundary conditions 
along the transverse directions, which are labeled by x and y coordinates, respectively. The 
space-time metric can be written as 

ds'' = -dt^ + Rlidxf + Rl{d^\ (1) 

and it will be convenient to define the volumes of the Neumann and Dirichlet directions 

Vj,^{Rj,Y^, VD^iRDY". (2) 

We take all dimensions to be larger than string scale 

Rn > 1, Rd> 1, (3) 

since one can apply a T-duality transformation to a small direction to make it large (we set 

a' = 1). 

An open string in such a background is labeled by liTv momentum numbers, rf^) winding 
numbers, the oscillator level and the two Chan-Paton factors identifying the D-branes on 
which the end points of the string are attached. Note that, unhke closed strings, there are 
no momentum modes along Dirichlet and no winding modes along Neumann directions. In 
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FIG. 1: Parallel D-branes in a compact space and open strings attached to them. 



the micro-canonical ensemble and for e ^ 1, the single string density of states is given by 
(see e.g. [35]) 



where the sum is over the integers < / < Rjj, Ph = 2\/27r is the Hagedorn temperature 
for open strings, eg = Rj), gi is a degeneracy factor^ which is approximately given by 
gi ~ 2Vol(5'^^-^)Z'^^-i for Z > 1 and 

Recall that n is the number of D-branes in the system and factor in (4) arises from the 
combinatorics of the Chan-Paton factors. 

The form of the single string density of states depends on the magnitude of the string 
energy e relative to the moduli parameter eo [35]: 



For e ^ eo only the first term contributes in the sum in (4) and in this limit the strings are 
energetic enough to wind around the Dirichlet directions. For example, this is the regime of 

interest for Rd = 0{1) and e ^ 1. In the opposite limit e <^ eo, the sum can be replaced by 
an integral and in this case the available energy is not large enough to excite the winding 
modes. 

The total density of states il(-E'), which counts the degeneracy with fixed energy E and 
unconstrained number of particles, is given by 



(4) 




(6) 





Indeed, gi counts the number of integer d£)-tuples whose length is equal to /, see [35]. 
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Using the integral representation of the delta function, Q{E) can be expressed as 



1 



(8) 



where 



Jo 



(9) 



To get the leading order asymptotic contribution one can use (4) in (9) to obtain 



(10) 



ia + E(3i 



where the temperature (3i is defined as 



(11) 



Note that F{a) is a regular function which vanishes as a — > ±00 and therefore the integral 
(8) is convergent. 

To calculate (8), we use a method which is similar to the complex temperature formalism 
[21-23]. Since F{a) is an analytic function in the complex a-plane (note that the points on 
the imaginary axis a = iEj3i are regular) one can deform the integral contour near the origin 
by extending it through the imaginary axis to circle around the point iE^^, as shown in 
figure 2. The analytic function F{a) can be written as the sum of two non-analytic functions 
F{a) — Fi{a) + ^2(0;), where Fi and F2 are the pieces containing the exponential and (— 1) 
in the denominator of (10), respectively. By expanding exp(Fi) in powers of Fi one can 
see that only the first term in the expansion contributes to the integral along the deformed 
contour. This is because, in evaluating all but the first term the contour can be closed from 
above by excluding the poles of Fi and F2 and the integral along the upper semi-circle at 
infinity vanishes. Therefore we have 



where the contour C is pictured in figure 2. This line integral can now be calculated by 
closing the contour from below; since there is no contribution coming from the lower semi- 
circle the integral equals to the sum of the residues at o; = iEf^i which we denote by fli{E): 




(12) 



Q{E) = Y,^i{E). 



(13) 
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FIG. 2: The deformed integration contour C, which is used to calculate the integral (8). 

It is easy to see that ^i{E) ~ exp(/5/-E'), therefore (13) offers an expansion scheme which can 
be seen to be equivalent to the complex temperature formalism [21-23]. In the following we 
concentrate on the calculation of VLq{E), the first term in this expansion . Since the number 
of poles depends on the radius of the Dirichlet directions, one should consider two different 
regimes depending on the magnitude of Rd- 

A. Small radius regime 

In the small radius regime, i.e. when Ro = 0{1), it is enough to consider only / = term 
in the sum in (12), and thus we have^ 



Expanding the first exponential and evaluating the residue at a = iE(3h term by term one 
finds [35] 



where Ji is the modified Bessel function of first kind. Compared to the closed string density 
of states, the main difference is the presence of the squareroot term in the exponential in 




(14) 




(15) 



^ As a convention, we define <j> to denote the line integral over a small circle around the singularity of the 
integrand, which is suitable for the calculation of the residue. 
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(15). Note that (15) is valid even when Rn ^ 1. Indeed for larger Rn the asymptotic 
expansion becomes more accurate. 

Although the residue can exactly be calculated in this case, one can also apply a saddle 
point approximation to extract the asymptotic behavior. After defining a = iE(3h + z y/fE 
the integral (14) becomes 

n{E) - e^'^^J^ (f dz e-*v^(--iA). (16) 



E 

For fE » 1, which is automatically satisfied in the small radius regime since Vd = 0{1), 
one can use a saddle point approximation about the point z = i which precisely gives (15). 



B. Large radius regime 

If Rd ^ 1, the I 7^ terms should also be taken into account in the sum in (12) while 
evaluating flo{E). This makes the exact calculation of the residue very difficult. However, 
under suitable conditions, a saddle point approximation can be used to get the asymptotic 
behavior. In terms of a new complex variable z — —{ia + PhE)/ E, the density function 

VLo{E) becomes 

no{E)^^ jdzexp{Ez + f/z + a{z)), (17) 
where a{z) is defined by 

^(-)-/E— |7- = E^^-'- (18) 

1^0 z + L /eo fc>o 

Note that a{z) is an analytic function near z — Oso that it can be expanded in a Taylor series 
as above. Using the fact that gi ~ 2Vol(5''^°~^)Z''^~^ for Z » 1 the expansion coefficients can 
be found as 

ak = {-l)''akfVD, (19) 

where Uk is a positive constant of order unity. With the help of this expansion and further 
defining 

K^^Jf/{E + a,), X^^f{E + a,), z^Ky, (20) 
(17) can be converted into 



X 



y + l/y+Y^ aky' 

k>2 



(21) 



where ak — ak/X. Assuming that A ^ 1, the saddle point approximation is apphcable to 
evaluate (21). Moreover, if cr^ <^ 1 for A; > 2, y = —i becomes the saddle point up to small 
controllable corrections. Under these assumptions the saddle point calculation gives [35] 

^^o(^) - _ ijv^)3/4 {^f^E + aofVn + 2^ f{E - a, fVo)) , (22) 

which is valid when A ^ 1 and a"jt <^ 1 for /c > 2. One can see that these two conditions are 
equivalent to 

E:^fVD + j. (23) 

For / > 1, (23) is identical to E > fVo- 

The contribution of the Tth pole fli{E) can be carried out in a similar fashion (^- variable 
should now be defined as —Ez — ia + (3iE). One can see that if the saddle point approxi- 
mation is applicable to D,o{E), i.e. when (23) holds, it is also pertinent in the evaluation of 
ili{E). A straightforward calculation then gives 

^ r M3/4 + aofVn + 2^ gi f {E - d, f Vn)) , (24) 

[gi{E - aJVon ^ ^ 

where do and di are two other constants of order unity. If one demands that Qq ^ ^i, an 

additional condition on E should be imposed 

E:^fRlVD or fR%, (25) 

which is stronger than (23) for / > 1. For do > 2, the first term becomes larger than the 
second term in the right hand side of the above inequality. As discussed in [35] , if (25) is 
not satisfied one has to change the form of the single string density of states u;{e) in the 
calculation of fl{E) according to (6). 

III. THE DILATON-GRAVITY BACKGROUND 

In this section we use dilaton-gravity equations sourced by open strings attached to D- 
branes to determine the string frame fields. Before deriving the field equations, let us discuss 
why such a setup might be relevant in the context of string gas cosmology. The original 
scenario proposed in [1] has been developed in [36] to include all higher dimensional excita- 
tions in string theory. As it is natural in such settings, all degrees of freedom are assumed 
in thermal equilibrium in a hot and dense state. Then, by applying the BV mechanism to 
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Dp-branes of different dimensions, it is argued in [36] that all higher dimensional branes 
with p > 3 annihilate in 10-dimensions and the remaining branes form a hierarchy of sizes 
of compact dimensions. Namely, 2-branes only permit a 5-dimensional subspace to grow, in 
which only 3-dimensions are allowed to expand by strings. 

As discussed in [36], even though BV mechanism works perfectly for higher dimensional 
branes, causality requires at least one brane per Hubble volume remaining. Although, a 
subsequent loitering phases can lead to a total annihilation [37], it is not very unnatural 
to assume that some D-branes survive even BV mechanism functions well. Besides, using 
Boltzmann equations the annihilation of D-branes in an expanding universe has been studied 
in [32] , which indicates that the B V mechanism may not work as efficient as one may think 
and higher dimensional D-branes may also fall out of thermal equilibrium surviving the 
annihilation. Therefore, it is worth to consider models with some D-branes left over in the 
system. 

At this point, one should notice that there is no conflict in assuming out of equilibrium 
D-branes and open strings in thermal equilibrium. In saying D-branes fall out of equilib- 
rium, one refers to D-brane anti-D-brane annihilation process, which happens only when 
branes physically intersect each other. For instance, if there exists only a single D-brane in 
the universe as it is usually assumed in brane-world models, then BV mechanism is inappli- 
cable, yet one can still think about open strings attached to that D-brane being in thermal 
equilibrium as it is studied in [34] or [35].^ 

As pointed out in the previous section if the net RR-charge of the universe is zero, then 
there must exist an equal number of D-branes and anti- D-branes. One then wonders whether 
the results of the previous section alter in such a modification. In [34], the thermodynamics 
of the system has been studied in this more general setup using the Boltzmann equations 
and (under suitable conditions) the number of states accessible with energy E is found to 
be the same (see eq. (19) in [34]). On the other hand, in a realistic situation one expects to 
encounter generic intersections rather than just parallel D-branes. Although it is possible 
to motivate the sole existence of parallel D-branes in the context of brane-world models, we 

^ It is interesting to note that when all directions are compact and at string size, the assumption of homo- 
geneity can be justified even in the presence of a single D-brane, since in thermal equilibrium open strings 
exist in a long string phase and they traverse the entire space many number of times before ending on the 
D-brane, which ensures homogeneity. 
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show at the end of this section that the main conclusion of this paper does not change if 
one considers generic intersections. 

Let us now start discussing the dilaton gravity equations. At weak string couphng gg — 
e*^ <^ 1, the action for the dilaton and the metric can be written as 



J d'^x^/^ [ e-^^{R + 4V0') + e-'t'Co 



(26) 



where Co is the effective Lagrangian for open strings and the dilaton dependence is explicitly 
singled out in (26). The field equations following from this action can be found as 



R^, + 2V^V,(/. -\[R + 4VV - 4(V0)^] g^, = e'^' T^,, 
R + 4V V - 4(V0)2 + ^e"^ Co = 0, 



(27) 
(28) 



where T^j^ is the energy momentum tensor 



T = 



1 



d 



•Co). 



(29) 



To proceed further one needs information about Co- It is noticed in [7] that imposing 
the conservation formula V^T^*^ = as the matter field equations, the contracted Bianchi 
identity can fix Co in terms of the energy-momentum content, at least in a cosmological 
setting. Taking 

= 0(i), gti^O, Tti^O, (30) 
where i labels a spatial direction, the consistency of (27) and (28) implies 



(31) 



which is precisely the Lagrangian for hydrodynamical matter. Therefore, given any con- 
served energy-momentum tensor, equations (27), (28) together with (31) determine the 
dynamics. 

In our case, the supposed D-brane configuration suggests 



ds' = -dt'' + e^^^'\dxf + e2^W(c^^)2, 



where the radii of the Neumann and Dirichlet directions are given as 



(32) 



i?jV — C ) Rd — 6 



C 



(33) 
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Assuming = and T^^, — disig{p,pN,PD), the equations (27) and (28) can be shown to 
imply 

+ i^0 = e^(^PA^ + ^PD-^p), (34) 

where dot denotes derivative with respect to t and 

K = (InB + (IdC - 20. (35) 

Compared to evolution equations obtained for closed strings that have no dilaton coupling, 
the system (34) has a very crucial difference: there appears the energy density p in the right 
hand side of B and C equations. Therefore, the "force" along a direction, which determines 
its cosmic evolution, is not only given by the pressure as in the case of closed strings but it 
is equal to sum of the pressure with energy density. 

The energy-momentum tensor for open strings on D-brane backgrounds can be derived 
from the entropy of the system determined in the previous section. Given S{E, Vn, Vd) one 
can find out the temperature and the pressures as 

^ -^^ P -TV P -TV (36) 

The densities that enter in the right hand side of (34) are given by 

E Pn Pd 

P^Y^ Piv = Pd = (37) 

where V = VnVd- Assuming that the cosmic evolution is adiabatic, i.e. dS = 0, implies 

E + dNBPN + doCPD = 0, (38) 

which is equivalent to conservation of energy-momentum tensor V^T'^'^ — 0. 
In the small radius regime, from (15) one finds 



Sc^PhE + 2JfE, (39) 
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and 




(40) 



which shows that the temperature is always smaller than the Hagedorn temperature and the 
pressures are very large. Not surprisingly, the pressures along Neumann and Dirichlet direc- 
tions turn out to be positive and negative, respectively, which is due to the presence/ absence 
of momentum and winding modes. 

Assuming that the universe starts out at the string radii Bq = Cq = Q with an expansion 
i?o, C'o > 0, one now has all the necessary information to determine the early cosmic evo- 
lution in this toy model. In principle one should specify 6 initial conditions for B, C and 
0, constrained by the last equation in (34). Also, the initial energy Eq (or constant entropy 
S) and the number of D-branes n in the system should be given. For £'o » 1, equations 
(39) and (40) can be used to determine the right hand side of (34) until say Rb ~ 3. After 
i^c > 3, one should consider the expressions in the large radius regime. 

We construct a simple analytic but approximate solution valid in the small radius regime 
as follows. Although the presence of D-branes breaks isotropy, we assume that initially 
Bq — Cq. We take n — 0{1) or so to avoid D-brane anti-D-brane annihilation (see the end 
of this section) and thus suppose 



Under these assumptions, one sees from (40) that T ~ to a very good approximation. 
Moreover, the pressure terms in the right hand side of (34) can be neglected compared to 
the energy density. Therefore, one effectively gets a pressure-less phase. Indeed, solving E 
from (39) 



one can see that it can approximately be treated as a constant (equal to S/ (3h) when / is not 
changing appreciably since S ^ f. This is consistent with an effective pressure-less phase. 

Ignoring the pressures, one can set B — C during the small radius regime. The system 
for B and can be rewritten in terms of the " conformal time" rj as 



Eq > ri^. 



(41) 




(42) 



2 




(43) 
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where 

e^^-'"^dr] = dt, (44) 

and prime denotes derivative with respect to rj. 

The first two equations give (j) — 3B + crj + (po, which can be used to get a single second 
order differential equation for B. There appears two more constants of integration in the 
solution of B one of which can be fixed in terms of the other constants using the constraint 
equation in (43). Choosing 77 = to be the initial time and setting also B{0) — we find 

B^X (e''^^ - 1) - ^ (45) 
= 3A (e^'^^ - 1) - 2c 77 + 00, (46) 



where 



= ^ e-"^° (47) 



and A is a free positive constant. By picking A to be positive, one makes sure that the 
proper time t increases with rj. For A > 1/9, the dimensions continuously expand where the 
expansion speed initially vanishes for A = 1/9. For A < 1/9, the initial expansion speed is 
negative giving a period of contraction which later turns into an ongoing expansion. These 
can be seen from the Hubble expansion parameter which can be determined as 

H = B = e-^^+^^B' = cexp [-3A(e3"'' - 1) - ct] + 20o] [sXe^""^ - l/s] . (48) 

Also note the initial Hubble rate 

//o = ce2^°[3A-l/3]. (49) 

Recall that this solution can be trusted during the small radius regime, i.e. until the time 
1]* such that e^^^*^ ~ 3. The change in dilaton during this period is also of order unity 
(f){rj*) — (f)(0) — 0{1). We give the plots of the functions exp{B), cj) and H ior X — 1/9 in 
figure 3. 

Having obtained the solution, one can now check the consistency of the background. 
Firstly, we should make sure that the prescribed initial conditions avoid Jeans instability, 
i.e. the gravitational collapse. For a black hole not to form in a region of size R in rf-spatial 
dimensions, the mass inside this region pW^ should yield a Schwarzschild radius"^ smaller 



Note that the Schwarzschild radius corresponding to a mass M in d-dimensions is given by = 
(GM)i/(rf-2). 
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FIG. 3: The plots of i? = e , H and for A = 1/9. We take c = 1 in the graphs. 



than R, which imphes 

< i-. (50) 

where G is the gravitational coupUng constant. Although for closed strings G = e'^'^, the 
coupling of open strings to gravity is determined hj G = e^. In our case since = 1 at 
t = 0, (50) gives 



which constraints the initial value of dilaton in terms of initial energy. Since the change in 
dilaton is of order unity, Jeans instability will be avoided at later times in the small radius 
regime, if it is avoided at r/ = with a margin. 

Secondly, we would like to check whether the assumption of thermal equilibrium can be 
justified in this model. This is an important constraint which is known to be problematic 
for closed strings [30]. Thermal equilibrium requires the interaction rate per string F to be 
larger than the expansion rate H: 



Although this condition has already been used in the literature in dilaton-gravity setting (see, 
e.g. [30, 31]), let us try to explain why it is still valid, since it is known that some standard 
results valid in Einstein gravity are modified in dilaton gravity. Equation (52) can be deduced 
from the Boltzmann equation for the number density n of a species, which roughly takes 
the following form in d-dimensions: h + dHn + Tn^ = . The Hubble term arises since the 
number density decreases with the scale factor as and the last term comes from the 

integration of the interaction cross sections. It is clear that in dilaton gravity the form of 
this equation is the same, consequently (52) can be used to justify thermal equilibrium. In a 
more fundamental setting, if one considers the Boltzmann equation for the number density 
in the phase space, then the Hubble term arises from the geodesic equation. In our case, the 




(51) 



T> H. 



(52) 
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conservation of the energy momentum tensor V^T'*'' = ensures that strings move on the 
geodesies in the string frame, which again shows that (52) can be used in dilaton gravity. 

We first evaluate F without paying attention to a possible "long" string dominance. In 
that case, the rate can be estimated as 

rshort--A^e'^ (53) 

where H is the number of strings per unit volume. To find Hi one can define the average 
number D{e) of strings carrying a fixed amount of energy e, which is given by [22] 

Since V = 0{1), the total number of open strings also equals to the number density, which 
can be found as 

C Die) de ~ \fjE, (55) 
Jo * 

where we used (6) and (15). Therefore, (53) gives 



short 



~ n 



v^e^"^, (56) 



where presently insignificant n dependence in the above formula is kept for future use. 

On the other hand, it is known that strings in a compact space near Hagedorn temper- 
atures can exhibit a long string phase where a single string can traverse the entire space 
several times. In that case, the interaction rate (53) should be modified to take into account 
this fact. Assuming a classical discrete model as discussed in [38, 39], the interaction rate 
for a single string should be proportional to its length times the total length of the rest of 
the strings 



(AA-l)L, (57) 



^ long ~ e 

where L is the average length of open strings. 

In our case, the average length L roughly equals to the average energy, which can then 
be found as 

E [E 

For n = C(l), the length becomes L ^ 1, which shows that the system is indeed dominated 
by long strings. Therefore the interaction rate can be estimated from (57) as 
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Although the interaction rate is greatly enhanced compared to short string estimate (56), 
equation (51) imphes that F <^ 1. 

Let us now check whether the assumption of thermal equilibrium can be justified. Using 
(49) and (59), the condition (52) at ?7 = implies that A should be fine tuned in a small 
neighborhood of 1/9. Recall that when A = 1/9 the initial expansion speed vanishes, so 
this is not surprising. Therefore, with a fine tuning it is possible to justify the assumption 
of thermal equilibrium initially. However, from (48) the Hubble parameter can be seen 
to increase and reach a maximum value at around crj ~ 0.4. The Hubble parameter then 
decreases a little bit, but at the end of the short radius regime, which corresponds to crj* ~ 1, 
it becomes roughly equal to its maximum value which can be determined from (48) as 

v^e"^"/'. (60) 

Looking at the solution, one finds that the dilaton is not changing appreciably during this 
interval. Therefore (52) implies 

^ <e'^°. (61) 



^2/3 

It is not possible to satisfy (61) along with (51) for E ^ 1 and e"^"^ ^ 1. Therefore, we see 
that even the initial conditions are fine tuned to avoid Jeans instability and justify thermal 
equilibrium, the open string gas most likely fall out of thermal equilibrium during the small 
radius regime. 

The only loophole in the above argument is that the estimated interaction rate V can 
actually be larger due to some numerical factors of vr's or sums over spin or momentum 
states etc. If these factors become large enough then it might be possible to avoid Jeans 
instability and assume thermal equilibrium for certain energies smaller than a critical value. 
However, this looks both difficult and unnatural. 

One may wonder whether this result might change if more generic intersecting configu- 
rations are considered instead of simply taking parallel D-branes. In that case, the spatial 
directions are divided into NN, DN, ND, or DD groups. As shown in [35], the single string 
density of states does not depend on ND and DN moduli, and the /-function in (5) becomes 

f = ^7^- (62) 

VDD 

Therefore, the total density of states also become independent of ND and DN moduli, and 
pressures along these directions exactly vanish. Furthermore, the pressures along NN and 
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DD directions are still ignorable compared to energy density, which shows that even in a 
generic situation involving intersecting branes one finds an effective pressure-less phase at 
string scale radii. Thus, the cosmic evolution and interaction rates do not change in the 
small radius regime and the assumption of thermal equilibrium is still questionable even 
intersecting D- branes exist. 

Another possible point of concern is the effect of closed strings on the evolution and the 
interaction rates. As shown in [34] , the Boltzmann equations imply that in equilibrium the 
average length and thus the average energy of closed strings is suppressed by the number 
of open strings jV, i.e. — Lo/M and Ec — Eg/M (see eq. (17) in [34]). Prom (55), we 
have M ^ 1 and thus the contributions of the closed strings on the cosmic evolution and 
interaction rates can safely be neglected due to the smallncss of their energy and the length. 

Although we have discovered that the background (45)-(46) is problematic, it is inter- 
esting to examine some further aspects of the solution. For instance, although the number 
of D-branes is taken to be small, one may want to inspect the issue of brane anti-brane 
annihilation. To be safe, one can simply require the Hubble expansion speed to be much 
larger than the peculiar brane speed, which can be assumed to be order one in string units. ^ 
Since the typical separation between D-branes is initially n~^^'^° , one thus demands 

l<//n-^/^°. (63) 

Using (60), this implies 



n 



2/dD 



^ <e^", (64) 

which can only be satisfied marginally along with the Jeans instability condition (64). 

As the Dirichlet directions expand, the expression for the entropy (39) looses its validity 
around Ru ~ 3 and one enters into a large radius regime. In the beginning of this regime 
we have^ Vyv — Vd and thus / > 1. From (22) the entropy is now given by 



S^(5hE + aofVD + 2Jf{E - ai f Vd), (65) 



^ Although strings are very energetic, most of their energies are stored in their masses and their motion 

can be viewed to be "non-relativistic" . D-branes embedded in such a bath of strings are not expected to 

move fast or stay motionless. 
^ Actually one expects Vn to be slightly larger than Vd, since (34) shows that the positive pressure along 

Neumann direction helps the expansion while the negative pressure along Dirichlet directions prevents it. 

Note that KB and KC terms play the role of a "velocity" dependent friction for K > 0. 



18 



where the condition (25) is assumed to be satisfied (note that (23) is automatically obeyed 
since / > 1). Since n — C(l), (25) is obeyed for E ^ 1. Prom (65), the temperature and 
pressures can be found as 



T \E-a,fVD' 



PM^TaofVD + T^f{E-a,fVD) - Ta^fvJ ^ ^ . , (66) 



PD = -T^Jf{E-a,fVD). 

As long as (25) is obeyed, we have E :$> f and E :§> /Vd, which shows that T ~ 1/Ph 
and E ^ Pjvj Pd- Thus, the pressures are still negligible in this regime and the solution 
(45)-(46) is valid until (25) is violated. Actually, since the dilaton is increasing in (46), there 
is also a possibility that the string coupling e'^ can become 0{1) and one enters in the strong 
couphng regime before even (25) looses its vahdity. Whether this happens or not depends 
on the initial conditions. In any case, as long as the equations (66) are trusted one ends up 
with an effectively pressureless evolution in this setup. 

At this point one may wonder why the arrow of time is chosen to yield an increasing 
dilaton and thus flow to a strong coupling regime. This choice actually dictated by the 
constraint equation in (34) as follows. According to the last equation, the sign of K cannot 
alter during the evolution. Initially the time fiow should be fixed to yield a positive K, since 
this choice produces "friction" terms in B, C and equations which avoids singularities. The 
initial conditions with K <Q give run-away solutions usually developing a naked singularity 
in a finite proper time and thus these are not suitable for the description of the universe 
following big-bang. In our case, the choice K > Q gives an increasing dilaton. 

One may also consider how the overall picture alters if the the number of D-branes n is 
large and comparable to initial energy n ~ Eq. In that case, (58) imphes that L — 0{1) and 
thus short strings dominate the system, which is due to presence of a large number of D- 
branes chopping up the long strings. Prom (40), one sees that the temperature is well below 
the Hagedorn temperature and pressures have the same order of magnitude with energy. 
From (40) defining an effective equation of state parameter u) as 

Pn^ujE, Pd^-ujE, ^ , , (67) 

one finds that initially a; ~ 0.1 if n ~ E'q. Since cu is changing in time it is difficult to solve 
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the field equations exactly or to utilize a simple approximation strategy. However, without 
even solving the system, one can observe that there is an obstruction which arises due to 
D-brane anti-D-brane annihilation. Equation (63) implies that to avoid brane anti-brane 
collisions the Hubble expansion speed must be large 

1 < Eq/'^" ~ n'/''^ < H. (68) 

On the other hand the interaction rate in the short string phase can be estimated from (56) 
as 

r ~ n e^"^ - 6^^° Eo^"^ < 1, (69) 

where the last inequahty follows from Jeans instability condition (51). Therefore, T <^ H 
and it seems difficult to avoid Jeans instability, assume thermal equilibrium and safely ignore 
D-brane annihilation process in this scenario. 

IV. CONCLUSIONS 

In this paper, wc consider a toy cosmological model dominated by open strings attached 
to D-branes and determine the corresponding dilaton-gravity solution. We use basic ther- 
modynamical properties of open strings on D-brane backgrounds to calculate the energy- 
momentum tensor, which turns out to be conserved under the assumption of adiabaticity. 
Consistent coupling of this conserved energy-momentum tensor to dilaton-gravity back- 
ground is determined using the contracted Bianchi identity. Contrary to closed strings, 
open strings attached to D-branes couple to dilaton and this alters the field equations in 
a significant way. Namely, the dynamical evolution of a direction is governed not only by 
the corresponding pressure as in the case of closed strings, but by the sum of pressure and 
energy density. Although the pressures are not small in string units, they can still be ne- 
glected compared to energy density and the early evolution is identical to a pressureless 
phase. Due to time depending dilaton, the metric functions differ from the usual matter 
dominated FRW cosmology. All directions and the dilaton tend to increase in the solution, 
which looses its validity until the energy density drops under a critical value or the strong 
coupling regime is reached. 

We check the self-consistency of the solution for a few issues. Firstly, avoiding Jeans 
instability imposes a constraint on the initial values of the dilaton and energy (51). Secondly, 
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to fulfill thermal equilibrium, initial conditions must be finely tuned such that the Hubble 
parameter becomes less than the interaction rate. However, even after this fine tuning, the 
Hubble parameter increases to a value (60) larger than the interaction rate, provided the 
condition (51) is assumed to avoid Jeans instability. Therefore, as for closed strings, the 
assumption of thermal equilibrium is questionable and one requires an understanding of non- 
equilibrium thermodynamics. On the other hand, the existence of suitable initial conditions, 
even though they are finely tuned, shows that the cosmology of open string gases may differ 
substantially from that of closed strings. 

In the context of string/brane gas cosmology, a complete and plausible scenario is still 
absent. Therefore, the role that can be played by D-branes and open strings attached 
to them is not clear in the big picture. It is somehow discouraging to observe that open 
strings suffer from a crucial puzzle encountered for closed strings. However, it looks like new 
features can arise that can modify the whole story. For instance, accepting the fine tuning 
in the initial conditions, the passage from thermal equilibrium to freezing out resembles the 
BV mechanism, so it would be interesting to dwell on this aspect by also including closed 
strings into the picture. Moreover, in such a scenario the Hubble parameter, which is initially 
vanishing, increases in time for a finite duration giving a decreasing Hubble radius. Such a 
behavior is desirable for the structure formation mechanism proposed in [40], therefore it is 
of interest to analyze revisions which can possibly arise by the inclusion of D-branes. 
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